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J£> ■ Abstract 

' We develop the formulation of turbulence in terms of the functional integral 

over the phase space configurations of the vortex cells. The phase space consists of 
Clebsch coordinates at the surface of the vortex cells plus the Lagrange coordinates 
of this surface plus the conformal metric. Using the Hamiltonian dynamics we find 
an invariant probability distribution which satisfies the Liouville equation. The 
violations of the time reversal invariance come from certain topological terms in 
effective energy of our Gibbs-like distribution. We study the topological aspects of 
^ ; the statistics and use the string theory methods to estimate intermittency. 
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1 Introduction 



The central problem of turbulence is to find the analog of the Gibbs distribution for the 
energy cascade. The mathematical formulation of this problem is amazingly simple. In 
inertial range we could neglect the viscosity and forcing and study the Euler dynamics of 
ideal incompressible fluid 

9 t V a = VpUJ a p - djl (1) 

Here 

U a/3 = d a Vf3 - d[3V a (2) 

is the vorticity field and 

h = P+\v 2 a (3) 
is the enthalpy which is eliminated from the incompressibility condition 

d a v a = (4) 

The key point is that the Euler dynamics can be regarded as a Hamiltonian flow in 
functional phase space. This geometric view was first proposed by Arnold in 1966 (see his 
famous book ||) and later developed by other mathematicians (see the references in the 



Moffatt's lecture in the 1992 Santa Barbara conference proceedings[10|). Here we derive 
the Hamiltonian dynamics from scratch using the conventional physical terminology. All 
the necessary computations are presented in Appendix A. Some of these results are new. 

The Hamiltonian here is just the fluid energy 



H = J d 3 r~v 2 a (5) 



and the phase space corresponds to all the velocity fields subject to the incompressibility 
constraint. The Poisson brackets between the components of velocity field 

WnWii)]=/AT,h-r)lVh-0^) (•> 

where 

Tap(r) = 5 aP 5(r)+d a d $ — (7) 

is the projection operator. The Euler equation can be written in a manifestly Hamiltonian 
form 

d t v a = f d 3 r'T a p (r - r') ^ 7 (r> 7 (r') = [v a , H] (8) 



The Liouville theorem of the phase space volume conservation applies here 

rt K(r),tf]=0 (9) 



Sv a (r) 

(Dv) = l[d 3 v(r) = const (10) 



r 
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The probability distributions P[v] compatible with this dynamics, must also be conserved, 
i.e. it must commute with the Hamiltonian 



[P[v},H} = (11) 

The Gibbs distribution 

P[v] = exp (-/?#) (12) 

is the only known general solution of the Liouville equation fllTI) . The above mentioned 



central problem of turbulence is to find another one. 

This formulation of turbulence is significantly different from the popular formulation 
based on the so called Wyld diagram technique ||. There, the functional integral is in 
place from the very beginning but the problem is to find its turbulent limit (zero viscosity). 
This functional integral involves time, so it describes the kinetic phenomena in addition 
to the steady state we are studying here. 

After trying for few years to do something with the Wyld approach I conclude that 
this is a dead end. The best bet here would be the renormalization group, which magi- 
cally works in statistical physics. Those critical phenomena were close to Gaussian, which 
allowed Wilson to develop the e expansion by rearranging the ordinary perturbation ex- 
pansion. 

There is no such luck in turbulence. The nonlinear effects are much stronger. The 
observed variety of vorticity structures with their long range interactions does not look like 
the block spins of critical phenomena. Moreover, there are notorious infrared divergencies, 
which make problematic the whole existence of the universal kinetics of turbulence. No! 
These old tricks are not going to work, we have to invent the new ones. 

If we are looking for something pure and simple this might be the steady state dis- 
tribution of vorticity structures. Here the geometric methods may allow us to go much 
further than the general methods of quantum field theory. Being regarded as a problem 
of fractal geometry rather than a nonlinear wave problem, turbulence may reveal some 
mathematical beauty to match the beauty of the Euler-Lagrange dynamics. 

This dream motivated the geometric approach to the vortex sheet dynamics in Q, 
where the attempt was made to simulate turbulence as the stochastic motion of the 
vortex sheet. This project ran out of computer resources, as it happened before to many 
other projects of that kind. However, the geometric tools developed in that work proved 
to be useful and we are going to use them here. 

Another false start: I tried to solve the Liouville equation variationally, with the 
Gaussian Anzatz with anomalous dimensions for the velocity field. The numbers came 
out too far from the experiment and it was hard to improve them. Similar attempts 



with the Gaussian Clebsch variables |TT[ also failed to produce the correct numbers. It 



became clear to me that velocity field fluctuates too much to be used as a basic variable 
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in turbulence.[] 

We encounter the same problem in QCD where the gauge field strongly fluctuates, 
and its correlations do not decrease with distance. The problem is not yet solved there, 
but some useful tools were developed. In particular, the Wilson loops (the averages of 
the ordered exponentials of the circulation of the gauge field) are known to be a better 
field variables. The Wilson loop is expected to be described by some kind of the string 
theory, though nobody managed to prove this. 

The dynamical equations for the Wilson loops as functionals of the shape of the loop 
were derived, and studied for many years [[J. The QCD loop equations proved to be very 
hard to solve, because of the singularities at self intersections. 

In my recent paper I derived similar equations for the averages of the exponentials 
of velocity circulation in (forced) Navier-Stokes equation. These equations have no singu- 
larities at self intersections, in addition they are linear, unlike the loop equations of QCD. 
This raised the hopes to find exact solution in terms of the string functional integrals. 

The theory developed below started as the solution the (Euler limit of the) loop 
equations. However, later I found how to derive it from the Liouville equation, whithout 
unjustified assumptions of the loop calculus. This is how I am presenting this theory here. 

2 Loop functional 

It is generally believed that vorticity is more appropriate than velocity for the description 
of turbulence. Vorticity is invariant with respect to Galilean transformations which shift 
the space independent part of velocity 

v a (r) v a (r - ut) + u a ; u a/3 (r) u a p(r - ut) (13) 

The correlation functions of velocity field involve the infrared divergencies coming from 
this part. Say, in the two- point correlation function this would be the energy density 

(va(r)v a (r')) = ™ - i<( Wa (r) - v a (r')) 2 ) (14) 
which formally diverges as 

H t°° 5 2 

— oc / dkk-i ~ Ls (15) 

V Jl/L 

according to Kolmogorov scaling |l[]. 

The infrared divergencies are absent in vorticity correlations. The corresponding 
Fourier integral is ultraviolet divergent due to extra factor of k 2 , but this is healthy. The 

1 This does not mean that the simple Gaussian models in velocity or Clebsch variables could not explain 
observed turbulence in finite systems. We are talking about idealized problem of isotropic homogeneous 
turbulence with infinite Reynolds number. 
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physical observables involve the vorticity correlations at split points, where the integrals 
converge. 

The complete set of such observables is generated by velocity circulations for various 
loops in the fluid 

T c [v] = j> dr a v a (r) (16) 

The constant part of velocity drops here after the integration over the closed loop, thus 
the circulation is Galilean invariant. The loop gets translated 

C^C-ut (17) 

but all the equal time correlations of the circulation stay invariant in virtue of translation 
symmetry. 

One can express the circulation in terms of vorticity via the Stokes theorem 



where S is an arbitrary surface bounded by C. In particular, for infinitesimal loop we 
would get the local vorticity. The moments of the circulation 

<(r c M) B ) = I [Dv)P[v] (T c [v)) n (19) 



all converge in the ultraviolet as well as in the infrared domain. The infrared convergence 
is guaranteed since these are surface integrals of vorticity, and the ultraviolet one is 
guaranteed since these are line integrals of velocity. 

This nice property suggests to study the distribution of the velocity circulation 

P C [T}= J[Dv]P[v}6(T-T c [v}) (20) 
It is more convenient to study the Fourier transform 

/oo r 
dTexp ( l7 r) P C (T) = / [Dv]P[v] exp (^cM) (21) 
-oo J 

We expect these functionals to exist in the turbulent limit unlike the distribution P[v}. 

The dynamical equation for these functionals (the loop equation) was derived in my 
previous work ||12|| . The time derivative of the circulation reads 

dtT c [v)= I dr a u aP (r)v p (r) (22) 
J c 

All the nonlocal terms reduced to the closed loop integrals of the total derivatives and van- 
ished. Being averaged with appropriate measure, the remaining terms in time derivative 
must vanish according to the Liouville equation 

{d t T c [v]e^(igTc[v})) = Q (23) 
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Note that this is not the Kelvin theorem of conservation of the circulation. The 
circulation is conserved in a Lagrange sense, at purely kinematical level 

jT c [v] = d t Y c [v\+jde 5 ^-v p {C{6)) = (24) 



In other words, the Euler derivative (p2|) is exactly compensated by the shift of every 
point at the loop by local velocity Now, according to the Liouville equation, each of 
these equal terms must vanish in average. 

The formal derivation goes as follows 

[Dv]P[v] exp (zjT c [v]) [t-yT c [v], H] (25) 

= J[Dv]P[v] [exp (z 7 r H),H] 

= - J [Dv] [P[v),B]exp(tjT c [v)) = 

and the physics is obvious: the average of any time derivative with time independent 
weight must vanish. 

Let us interpret in these terms the solution for ^c{l) found in the previous paper 

*c(7) = = <f dr a Tp (26) 

In virtue of linearity of the Liouville equation it suffices to check the Fourier transform 

exp(z^R a pE a p) (27) 
In terms of the velocity field this corresponds to global rotation 

v a = R a pr p; to a p = 2R a p (28) 
The corresponding Poisson brackets reads 

[T c [v],H] = i dr a u a pvp = AT Jai R a pRp 1 = (29) 
J c 

The sum over tensor indexes vanished by symmetry. 

This solution is always present in fluid mechanics due to the conservation of the angular 
momentum. Unfortunately, it has nothing to do with turbulence, contrary to the hopes 
expressed in my previous work. 

Let us also note, that the Gibbs solution does not apply here. One could formally 
compute the loop functional for the Gibbs distribution, but the result is singular^] 



*c(7) = ex P (y~Jpf c dr <* f c dr '^ 3 



(r - r') (30) 



2 This is yet another advantage of the loop functional: the Gibbs solution for the loop functional does 
not exist, which forces us to look for alternative invariant distributions. 
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Should we cut off the wavevectors at k ~ — this would become 

ro 

* c ( 7 )«exp^--^jfjdr|) (31) 

This is so called perimeter law, characteristic to the local vector fields. Clearly this is 
not the case in turbulence, as velocity field is highly nonlocal. Also, the odd correlations 
of velocity, such as the triple correlation, which are present due to the time irreversibility, 
would, in general make the loop functional complex. 



3 Vortex dynamics 

Let us study the dynamics of the vortex structures from the Hamiltonian point of view. 
We shall assume that vorticity is not spread all over the space but rather occupies some 
fraction of it. It is concentrated in some number of cells D, L of various topology moving 
in their own velocity field. 

Clearly, this picture is an idealization. In the real fluid, with finite viscosity, there 
will always be some background vorticity between cells. In this case, the cells could be 
defined as the domains with vorticity above this background. The reason we are doing 
this is obvious: we would like to work with the Euler equation with its symmetries. 

We shall use two types of tensor and vector indexes. The Euler (fixed space) tensors 
will be denoted by Greek subscripts such as r a . The Lagrange tensors (moving with 
fluid) will be denoted by latin subscripts such as p a . The field X a (p) describes the instant 
position of the point with initial coordinates p a . The transformation matrix from the 
Lagrange to Euler frame is given by d a X a (p). 

The contribution of each cell D to the net velocity field can be written as follows 

^-^L^w^m (32) 

where 

q 7 (p) = n a {p)d a x 1 { P ) (33) 

is the vorticity vector in the Euler frame. The vorticity vector Sl a (p) in the Lagrange 
frame is conserved 

d t n a (p) = (34) 

The physical vorticity tensor u) a p = d a V[3 — d[3V a inside the cell can be readily computed 
from velocity integral. The gradients produce the S function so that we get 

^ M)= K:3 "''^ M (35) 
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or in Euler frame, inverting p = X 1 (r) 

Wapir) = g e abc d x p a d^p b d u p c e a/37 Q f (p)d f X 1 = e abc d a p a (r)dpp b (r) fi c (p(r)) (36) 

The inverse matrix d a p a = (9 a X a ) _1 relates the Euler indexes to the Lagrange ones, as 
it should. These relations between the Euler and Lagrange vorticity are equivalent to the 
conservation of the vorticity 2-form 

n=J2 uj aP (x)dx a Adx p = Y, n ab (p)d P a a d P b ; n ab = e abc n c (37) 

a<(3 a<b 

which is the Kelvin theorem of the conservation of velocity circulation around arbitrary 
fluid loop. 

The field X a (p) moves with the flow (the Helmholtz equation) 

d t X a (p)=v a (X(p)) (38) 

In Appendix B we derive the Euler equation from the Helmholtz equation and study the 
general properties of the former. We show that this is equivalent to the Hamiltonian 
dynamics with the (degenerate) Poisson brackets 

[X a (p),X p (p')} = -^(P-P')e^^ (39) 

The Hamiltonian is given by the same fluid energy, with the velocity understood as func- 
tional of X(.) The degenaracy of the Poisson brackets reflects the fact that there are only 
two independent degrees of freedom at each point. This leads to the gauge symmetry 
which we discuss below. 



The Hamiltonian variation reads 

5H 



5X a (p) 



e afiy v p (X(p))n y (p) (40) 



This variation is orthogonal to velocity, which provides the energy conservation. It is also 
orthogonal to vorticity vector which leads to the gauge invariance. The gauge transfor- 
mations 

5X a (p) = f(p)Q a (p) (41) 
leave the Hamiltonian invariant. These transformations reparametrize the coordinates 

5p a = f( P )n a (p) (42) 

The vorticity density transforms as follows 

5Q a = -Q a d b {fQ b ) + fQ b d b Q a + Q b d b {fQ a ) = 2ftt b d b Q a (43) 



8 



The first term comes from the volume transformation, the second one - from the argument 
transformation and the third one - from the vector index transformation of Q so that 



d 3 pQ a d a = inv (44) 

The identity 

d a tt a = (45) 

was taken into account. We observe that these gauge transformations leave invariant the 
whole velocity field, not just the Hamiltonian. 

The vorticity 2-form simplifies in the Clebsch variables 

fT = e abc d b Mp)dcMp) (46) 
Q = d(pi A d(j) 2 = inv (47) 

The Clebsch variables provide the bridge between the Lagrange and the Euler dynamics. 
By construction they are conserved, as they parametrize the conserved vorticity. The 
Euler Clebsch fields $i(r) can be introduced by solving the equation r = X(p) 

Mr) = &(x- 1 (r)) (48) 

However, unlike the vorticity field, the Clebsch variables cannot be defined globally in 
the whole space. The inverse map p = X~ l (r) is defined separately for each cell, therefore 
one cannot write v a = &id a § 2 + d a § 3 everywhere in space. Rather one should add the 
contributions from all cells to the Poisson integral, as we did before. 

This explains the notorious helicity paradox. The conserved helicity integral 

H = f d 3 pv c (p)Q c (p) (49) 

where 

V a {p) = (f>ld a fo + d a (p 3 ; <9 2 3 = -<9 a (01<9 a 2 ) (50) 

is the initial velocity field (to be distinguished from the physical velocity field v a (r) which 
cannot be paramatrized globally by the Clebsch variables). 

The helicity integral for the finite cell could be finite. It can be written in invariant 
terms of the vorticity forms 

Hv= I d z pe ahc d a <\) X d h fad c fo — ( d(f>i A d0 2 A d<p 3 = f 3 #1 A dcfo (51) 

JD JD JdD 

from which representation it is clear that it is conserved 

U v = I $3 d&x A cM> 2 (52) 

JdD 

Our gauge transformation leave the Clebsch field invariant 

5UP) = f{p)^ a {p)d a Up) = (53) 
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The velocity integral in Clebsch variables reduces to the 3-form 

dXy A d<fii A d(j)2 



v a (r) = -e a/37 dp J d 4?r|r _ x| (54) 

This gauge invariance is less than the full diffeomorphism group which involves arbi- 
trary function for each component of p. This is a subtle point. The field Q a (p) has no 
dynamics, it is conserved. However, the initial values of fl a can be defined only modulo 
diffeomorphisms, as the physical observables are parametric invariant. So, we could as well 
average over reparametrizations of these initial values, which would make the parametric 
invariance complete. 

Another subtlety. The equations of motion, which literally follow from above Poisson 
brackets describe the motion in the direction orthogonal to the gauge transformations, 
namely 

dt X a (p) = - " a ffify (p) ) v p (X(p)) = v a (X(p)) + f(p)n a (p) (55) 

The difference is unobservable, due to gauge invariance. We could have defined the 
Helmholtz equation this way from the very beginning. The conventional Helmholtz dy- 
namics represents so called generalized Hamiltonian dynamics, which cannot be described 
in terms of the Poisson brackets. The formula ( |I0|) cannot be solved for the velocity field 
because the matrix 

ttap{p) = ea/3 7 ^ 7 (p) (56) 

cannot be inverted (there is the zero mode fig(p)). The physical meaning is the gauge in- 
variance which allows us to perform the gauge transformations in addition to the Lagrange 
motion of the fluid element. 

Formally, the inversion of the Q— matrix can be performed in a subspace which is 
orthogonal to the zero mode. The inverse matrix f2^ 7 in this subspace satisfies the equation 

n^n* = 5 ai - %^ (57) 

which has the unique solution 

^ 7 = ~e a ^ (58) 

This solution leads to our Poisson brackets. 

Now we see how the correct number of degrees of freedom is restored. The Hamilto- 
nian vortex dynamics locally has only two degrees of freedom, those orthogonal to the 
gauge transformations. We could have obtained the same Poisson brackets by canonical 
transformations in the vorticity form from the Clebsch variables to the X— variables. This 
canonical transformation describes the surface in the X— space. The vorticity form at this 
surface can be treated as a degenerate form in the 3-dimensional space. 
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One may readily check the conservation of the volume element of the cell 

d (Xi,X 2 ,X 3 ) d(Xi,X 2 ,X 3 ) 
d t r = ^T7 \~ d aV a {X) = (59) 

0{P1,P2,P3) d{pi,p 2 ,p3) 

Using the formulas for the variations of the velocity field derived in Appendix B one 
may prove the stronger statement of the phase volume element conservation (Liouville 
theorem) 

(DX) = Y[d 3 X(p) = const (61) 



4 Vortex statistics 



Let us recall the foundation of statistical mechanics. The Gibbs distribution can be for- 
mally derived from the Liouville equation plus extra requirement of multiplicativity. In 
general, any additive conserved functional E could serve as energy in the Gibbs distribu- 
tion. 

The physical mechanism is the energy exchange between the small subsystem under 
study and the rest of the system (thermostat). The conditional probability for a subsystem 
is obtained from the microcanonical distribution dT'dT5(E'+E~E ) for the whole system 
by integrating out the configurations dV of the thermostat. 

The corresponding phase space volume e s ^ E '^ = J dT'S(E' + E — E ) of the thermostat 
depends upon its energy E' = E — E where the contribution — E from the subsystem 
represents the small correction. Expanding S(Eq — E) = S(Eq) — [3E we arrive at the 
Gibbs distribution. 

In case of the vortex statistics we may try the same line of arguments. The important 
addition to the general Gibbs statistics is the parametric invariance. Reparametrizations, 
or gauge transformations, are part of the dynamics, as have seen above.0 So, the Gibbs 
distribution should be both gauge invariant and conserved. 

The net volume of the set of vortex cells V = J2iV{Di) where 

V(D) = I d 3 p 9 ( XuX2 ' X ^ — [ dX\ A dX 2 A dX 3 (62) 

JD 0(j>i,p2,p3) JD 

is the simplest term in effective energy of the Gibbs distribution. It is gauge invariant, 
additive and positive definite.lt is bounded from above by the volume of the system, so 
that there could be no infrared divergencies.^ 

3 This makes so hard the numerical simulation of the Lagrange motion. The significant part of notorious 
instability of the Lagrange dynamics is the rcparamctrization of the volume inside the vortex cell. 

4 The Hamiltonian does not exist in the turbulent flow because the energy spectrum diverges at small 
wavevectors. The net Hamiltonian grows faster than the volume of the system which is unacceptable for 
the Gibbs distribution. 
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The mechanism leading to the thermal equilibrium is quite transparent here. The 
volume of a little cell surrounded by the large amount of other cells, would fluctuate due 
to exchange with neighbors.QThese are the viscous effects, in the same way as the energy 
exchange mechanisms in the ideal gas were the effects of interaction. The relaxation time 
is inversely proportional to the strength of interaction (viscosity in our case). 

We must take these effects into account in kinetics, but the resulting statistical distri- 
bution involves only the energy of the ideal system. This was the most impressive part 
of the achievement of Gibbs and Boltzman. They found the shortcut from mechanics to 
statistics, avoiding the kinetics. All the interactions are hidden in the temperature and 
chemical potentials. 

Mechanically, the cells avoid each other as well as themselves and preserve their topol- 
ogy but the implicit viscous interactions would lead to fluctuations. Even if we start from 
one spherical cell it would inevitably touch itself in course of the time evolution. At the 
vicinity of the touching point the viscous effects show up, which break the topological 
conservation laws of the Euler dynamics. The result could be a handle, or the splitting 
into two cells. After long evolution we would end up with the ensemble of cells Di with 
various number hi of handles. 

The related subject is the vorticity vector field fl a (p) inside the cells. In the Euler dy- 
namics it is conserved, but the viscosity-generated interaction would lead to fluctuations. 
The invariant measure is 

(DO) = l[d 3 n(p)6[d a n a (p)} (63) 



In terms of the Clebsch variables (46) the measure is simply 

(D<P) = Y[d 2 <P(p) (64) 
p 

as these are canonical hamiltonian variables. As discussed above, there are no global 
Clebsch variables in Euler sense. These Clebsch variables are defined separately in each 
cell and the net velocity field is a sum of contributions from all cells rather than the single 
Clebsch-parametrized expression $ 1 5 a $ 2 + d a ^ 3 . 

What could be an effective energy for the vorticity? The helicity integral is excluded 
as a pseudoscalar, besides, it is nonlocal, like the Hamiltonian. We insist on parametric 
invariance and locality in a sense that the cell splitting and joining do not change this 
energy. The Clebsch variables are defined modulo additive constants, they could be 
multivalued in complex topology, therefore we have to use the vorticity field itself. The 
generic Q— invariant in d dimensions is 



J d d psfM Q ab (65) 



5 The volume, as well as any other local functional of the cell does not change at splitting/joining, there- 
fore these processes go with significant probability. For the nonlocal functionals, such as the hamiltonian, 
there are long range interactions, which makes the exchange process less probable. 
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In even dimensions d = 2k this invariant reduces to the Pfaffian 

f d 2k p^fdet n ab = j ft A ft ... A ft = j ^ (66) 
In two dimensions it would be simply the net vorticity of the cell 

d 2 P Q 12 (67) 



Id ' 



However, in odd dimensions it vanishes so that there is no ft— invariant. We see that there 
is a significant difference in the vortex statistics in even and odd dimensions. 

Another interesting comment. For odd k — | the ft— invariant is an odd functional 
of ft which breaks time reversal invariance. The vorticity stays invariant under space 
reflection but changes sign at time reversal. This simple local mechanism of the time 

irreversibility is present only at d — 4k + 2,k — 0,1, In three dimensions we live in it 

is absent. 

Let us turn to the boundary terms. The boundary of the cell S = dD is described by 
certain parametric equation 

S: Pa = i2«(£i,6) (68) 

Clearly, in our case this is a self- avoiding surface. Here we could add the following local 
surface terms to the energy 

E<t> = E / d2 Z ( a V9 + by/ggvdjfadjfa) (69) 

where 

9iA0 = dMOdjRaiO (70) 

is the induced metric. 

There is also a topological term for each non-contractible loop L of dD 

E e = Q J fadfa (71) 

This is simply the velocity circulation around such loop. Assuming continuity (i.e. van- 
ishing) of vorticity at the boundary of the cell, this circulation can be also written as the 
integral in the external space. Then is it obvious that this integral does not depend upon 
the shape of the loop L, it is given by the invariant vorticity flux through the cross section 
E of the corresponding handle 



E e = e 



J^VL; L = <9£ (72) 



This term breaks the time reversal! This is the only possible source of the irreversibility 
in this theory in three dimensions. 
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The following observation leads to crucial simplifications. The only term which de- 
pends upon the Lagrange field X a (p) is the volume term, which in fact is the functional 
of the bounding surface X(dD) 

V= f dXi A dX 2 A dX 3 = [ X 3 dX 1 AdX 2 (73) 

JD JdD 

The rest of the terms in the effective energy also depending only upon the boundary, 
this lowers the dimension of our effective field theory. We are dealing with the theory 
of self-avoiding random surfaces rather that the 3-d Lagrange dynamics. With some 
modifications the methods of the string theory can be applied to this problem. 

The invariant distance in the X(£) functional space is 

| \5X\ | 2 = f gD <PZy/g(6X a {Q) 2 (74) 

where g as usual stands for the determinant of the metric tensor. One may check that 
the corresponding volume element 

(DX) =U6X (75) 
ZedD 

is conserved in the Euler-Helmholtz dynamics as well as the complete volume element. 
The key point in this extension of the Liouville theorem is the observation that the matrix 
trace of the functional derivative vanishes 

for arbitrary p, p' , including the boundary points. The metric tensor g^ does not introduce 
any complications, as it is X independent. 

This metric is the motion invariant in the vortex dynamics. In the statistics, according 
to the general philosophy these invariants become variables. We see that the field R a {0 
enter only via the induced metric, which allows us to introduce the latter as a collective 
field variable. 

One has to introduce the functional space of all metric tensors with the Polyakov 
distance 

I \Sg\ | 2 = / d 2 i^~g (5g tl 5g kl ) Ug» g kl + Bg*g* 1 ) (77) 
J s 

The parametric invariance can be most conveniently fixed by the conformal gauge 

9d0 = 9ij(0e aip(0 (78) 

where g is some background metric parametrizing the surface with given topology. Unlike 
the internal metric, the background metric does not fluctuate. We have the freedom to 
choose any parametrization of the background metric. 
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The effective energy which emerges after all computations of the functional jacobians 
associated with the gauge fixing readsf] 

E v = ^ / d 2 ijg - QR^p + iiz av ) (79) 

where R is the scalar curvature in the background metric. 

The parameters Q, a should be found from the self-consistency requirements. In case 
of the ordinary string theory in d dimensional space the requirement of cancellation of 
conformal anomalies yields 0] 



VT=d-V25^d 25-d 

a= — — ;0 = V~ (80) 

In three dimensions a is a complex number, which is fatal for the string theory. Fortu- 
nately this formula does not apply to turbulence, because the dynamics of the X field is 
completely different here. Later we speculate about the values of these parameters. 

To summarize, the total phase space volume element of our string theory 

dT= J] (DX)(D<j>){D<p) (81) 
cells 

and the total effective energy 

E= lv + E 4 , + E v + £ E @ \ (82) 
cells \ loops 



The grand partition function 

N 



Z = Y exp (-fiN - A Y *k\ f dT exp {-(3E) 



33) 



The interaction between cells comes from the excluded volume effect. 
The vorticity correlations are generated by the loop functional [|12f 

* c (7) = <exp(«rcrM)) (84) 

where 

r c [v] = f dr a v a {r) (85) 
J c 

is the velocity circulation. Our solution for the loop functional reads 

*o( 7 ) = Z- 1 Y exp (-fiN -\Y,hi) f dT ex P (-0E + % jT c [v}) (86) 

N,hi \ i=l / 



6 We cut some angles here. This effective measure was obtained [^| with some extra locality assump- 
tions in addition to the honest calculations. These assumptions were never rigorously proven, but they 
are known to work. All the results which were obtained with this measure coincide with those obtained 
by the mathematically justified method of dynamical triangulations and matrix models. 
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where / implies that the cells also avoid the loop C. 

There are now various topological sectors. In the trivial sector, the loop can be 
contracted to a point without crossing the cells. Clearly, circulation vanishes in this 
sector. In the nontrivial sectors, there is one or more handles encircled by the loop, so 
that the circulation is finite. Here is an example of such topology 



The circulation can be reduced to the Stokes integral of the vorticity 2-form over the 
surface Sc encircled by the loop C. Only the parts Sj = Sc H A passing through the 
cells contribute. In each such part the vorticity 2-form can be transformed to the Clebsch 
coordinates which reduces it to the sum of the topological terms 



In presence of the same terms in the effective energy the loop functional would be com- 
plex, as the positive and negative circulations would be weighted differently. This is 
manifestation of the time irreversibility. 

The interesting thing is that the interaction between the X field and the rest of 
field variables originates in the topological restrictions. The circulation reduces to the 
net flux from the handles encircled by the loop in physical space (this is a restriction 
on the X field). Also, the requirement that the loop is avoided by cells imposes another 
restriction on the X field. There are no explicit interaction terms in the energy. The whole 
dependence of the loop functional on the shape of the loop C comes from the excluded 
volume effect. The implications of this remarkable property are yet to be understood. 

Let us now check the Liouville equation. The time derivative of the circulation reads 




T C = E / fi = E / #1 A #2 = E / <M02 



(87) 



(see(HD) 




(88) 



But the vorticity vanishes at the loop, because the cells avoid it! Therefore 



(d t T c [v] exp (iyTc[v])) = 



(89) 



which is the Liouville equation. 
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For the readers of the previous paper |12| let us briefly discuss the loop equation. The 
essence of the loop equation is the representation of the vorticity as the area derivative 
acting on the loop functional 

(90) 



r) = 



7 Sa a p(r) 

The formal definition of the area derivative in terms of the ordinary functional derivatives 
was discussed before [|], [L2|. The geometric meaning is simple: add the little loop to 
the original one and find the term linear in the area 5a a p enclosed by this little loop. 
The parametric invariant functionals like this one can always be regularized so that the 
variation would start from 5a a p. The area derivatives of the length and the minimal area 
inside the loop were derived in j4|. 

The velocity operator is related to vorticity by the Poisson integral 



-dp I d'R 



An\r - R\ 



(91) 



The geometric meaning is as follows. The vorticity operator adds the little loop SCr^r at 
the point R off the original loop C. By adding the couple of straight line integrals 



W(r, R) = exp yl J 



W~\R,r) 



(92) 




we reduce this to one loop of the singular shape 



C — {C r>r , L r< R, SCr,R, Lr )T ) (93) 

Tc + = I> (94) 



which can be obtained from original one as certain variation||12||. Then the loop equation 
is simply 

' dr a u a p(r)vp(r)q c (~f) = (95) 



ic 

It was shown in [|EJ that these two operators commute as they should. 

In our case this equation is satisfied in a trivial way. Regardless the subtleties in the 
definition of u, v, as long as these operators act on the vortex sheet circulation as they 
should, they insert the vorticity and velocity at the loop. The rest of the argument is the 
same as in the Liouville equation. 

Note that this trick would not solve the Hopf equation for the velocity generating 
functional 

H[J] = (exp (| d 3 rJ a (r)v a (r^j) (96) 
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The Hopf equation would require 

( J d 3 r'T a(3 (r - r') LO^{r')v 7 {r') exp (J d 3 rJ a (r)v a (r)^j) = (97) 

which we cannot satisfy since the volume integral / d 3 r' would inevitably overlap with 
cells where vorticity is present. 

Apparently we found invariant probability distribution for vorticity but not for the 
velocity. The velocity distribution may not exist in the infinite system, due to the infrared 
divergencies. In practice this would mean that the velocity distribution would depend of 
the details of the large scale energy pumping, but vorticity distribution would be universal. 

The correlation functions of vorticity field can be obtained from the multiloop func- 
tional by contracting loops to points. In case when the loops encircle one or more handles 
there would be nonvanishing correlation function. Here is an example of the topology with 
the two point correlation which also has nontrivial helicity because of the knotted handle. 




5 Discussion 

So far we do not have much to tell to the engineers. Even if this statistics of the vor- 
tex structures will be confirmed by further study, we shall face the formidable task of 
computing the correlation functions. Let us speculate what could come out of these 
computations. 

The qualitative picture of intermittent distribution of vorticity will be the same as in 
the multifractal models flOU - In fact these models inspired our study to some extent. Our 
basic idea is that the part V of the volume occupied by vorticity fluctuate. In numerical 
and real experiments the high vorticity structures were clearly seen. The cells take 
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the shape of long sausages rather than spheres, which does not contradict our general 
philosophy but still lacks an explanation. 

Let us try to estimate the intermittency effects in the vortex cells statistics. Let us 
contract the loop to a point r around some handle. What we get in the limit can be 
expressed in terms of the vertex operator of the string theory 

Tc oc f d%^ f d 2 Vo ^6 3 (X(0 - r)S 3 (X( V ) - r'); r' - r (98) 

J S tJ s 

The points £o and rjQ are mapped to the same point r' — > r in physical space: this is 
the handle strangled by the loop. The important detail here is the factor ^fg = e aip 
corresponding to the metric tensor at the surface. 

The properties of such metric tensors were studied in the string theory H ||. The 
moments of Tq would behave as 

(r£) oc r A{n) ; A(n) = \a{na + Q) (99) 

where ro is the short-distance cutoff. The parameters a, Q are to be found from the 
selfconsistency conditions. 

In case of the turbulence theory we know that the third moment has no anomalous 
dimension, due to the Kolmogorov's | law. This implies that 

Q = -3a (100) 

after which we exactly reproduce the anomalous dimensions of so called Kolmogorov- 
Obukhov intermittency model. 

They assumed in 1961 the log-normal distribution of the energy dissipation rate as 
a modification of the Kolmogorov scaling. Various multifractal models generalizing this 
distribution, and the corresponding experimental and numerical data are discussed in 
proceedings of the 1991 conference in Santa Barbara ||10|| . 

The physics of the fractal dimensions in our theory is almost the same as in the 
Kolmogorov- Obukhov model. The local energy dissipation rate at the edge can be esti- 
mated as oc e 3a,fi . This quantity is, indeed, an exponential of the gaussian fluctuating 
variable. The variance is proportional to logr since it is the two dimensional field with 
the logarithmic propagator (how could they have guessed that!). 

Strictly speaking, the dimensions of powers of Tq were never measured. People con- 
sidered the moments of velocity differences instead. In principle, the potential part of v 
which is present in these moments may change the trajectory, so we must be cautious. 

The point is, the potential part has completely different origin in our theory. It comes 
from the nonlocal effects, involving all the scales, including the energy pumping scales. In 
short, this part is infrared divergent. The vorticity part which we compute, is ultraviolet 
divergent, it it determined by the small scale fluctuations of the vortex sheet, represented 
by the Liouville field. 
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In absence of direct evidence we may try to stretch the rules and estimate our inter- 
mittency exponents from the moments of velocity. I would expect this to be an upper 
estimate, as intermittency tends to decrease with the removal of the large scale effects. 
In my opinion, there is still no direct evidence for the anomalous dimensions in vorticity 
correlators. It would be most desirable to fill this gap in real or numerical experiments. 

Let us stress once again, that above speculations do not pretend to be a theory of 
turbulence. Still they may give us an idea how to build one. 

In statistical mechanics the Gibbs distribution was the beginning, not the end of the 
theory. If this approach is correct, which remains to be seen, then all the work is also 
ahead of us in the string theory of turbulence. The string theory methods should be fitted 
for this unusual case, and, perhaps, some scaling and area laws could be established. I 
appeal to my friends in the string community. Look at the turbulence, this is a beautiful 
example of the fractal geometry with extra advantage of being guaranteed to exist! 
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A Poisson brackets for the Euler equation 

One may derive the Poisson brackets for the Euler equation by comparing the right side 
of this equation with 

[v a (r),H] = J d"r'[v a (r),v^r')] J ^ = J d 3 r' [v a (r), v p {r')] v p {r') (101) 

Comparing this with the Euler equation we find 

[v a (r),vp(r')} = T au (r - r')cu^(r') + . . . (102) 

where dots stand for the dp terms which drop in the above integral. These terms should 
be restored in such a way that the Poisson brackets would become skew symmetric plus 
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they must be divergenceless in r' as well as in r 



' K(r), V/3 (r')] = ^K(r),^(r')]=0 (103) 



dr a ' dr^ 

The unique solution is given by the formula in the text. 

It is also worth noting that one could derive the Poisson brackets for velocity field 
using Clebsch variables, which are known to be an ordinary (p, q) pair. The velocity field 
is represented as an integral 

v a (r) = J d 3 r'T a(S (r - r')$ i(r , )<9 /3 $ 2 (r') (104) 

and the corresponding vorticity is local 

u a p(r) = e'id&dpQj (105) 

The Euler equations are equivalent to the following equations 

6 t + vpdp^ = (106) 

which simply state that the Clebsch fields are passively advected by the flow. These 
equations have an explicit Hamiltonian form 

r tt 

h = S Wj (107) 

where it is implied that the Hamiltonian is the same, with velocity expressed in Clebsch 
fields. 

Now the Poisson brackets for the velocity fields can be computed in a standard way 

which again yields the formula in the text. This derivation is not as general as the previous 
one, since there are some flows which cannot be globally described in Clebsch variables. 



B Helmholtz vs Euler dynamics 

Let us study the variation of the vortex cell velocity field. Our first objective would be 
to show that the time variation according to the Helmholtz equation 

5X a (p) = dtv a (X(p)) (109) 

reproduces the Euler equation 

8v a (r) = dt Mr)u; a/3 (r) - d a h{r)) (110) 
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Simple calculation of the variation of the initial definition with integration by parts 
yields 

6v a (r) = {e aPu d, - e aPl d v ) dp [ d 3 p ^(p)^r(p) (m) 

Jd 47r|r — A{p)\ 

Now we use the identity 

e a /3vd-y — e a/ 3 7 d u = e^ v d a — e aiu dp (112) 

(the difference between the left and the right sides represents the completely skew sym- 
metric tensor of the fourth rank in three dimensions, which must vanish). This gives us 
the following two terms in velocity variation 



5v a (r) = -d a 5h(r) + e aiu j d'pSX^Q^S^r - X(p)) (113) 

where 

6hir) = e ^J dp An\r-X(p)\ (114) 

The first term is purely potential, it can be reconstructed from the second term by 
solving d a Sv a (r) = 0. Now we see that the Helmholtz variation reproduces the Euler 
equation with the enthalpy 

h(r) = e^d p J d p 4?r|r _ x(p)| (H5) 

At the same time we see that the Helmholtz variation is defined modulo gauge trans- 
formation 

8X v {p) =s> 5X u (p) + f{p)n v {p) (116) 
which leaves the velocity variation invariant. 
Note that the functional derivative 

SXM ~ {ea " A " e ^ v) d "4n\r-X(p)\ (U7) 
is a traceless tensor in a, v . This is sufficient for the volume conservation 

which is the Liouville theorem for the vortex dynamics. 
We found the following Poisson bracket for the X field 

[X a (p),X p (p')} = -S"(p-p')e aPl ^ (119) 

The equations of motion corresponding to these Poisson brackets read 

dt X a (p) = [X a (p),H] = J d 3 p> [X a {p),X p {p')\ -^-^ (120) 
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Let us compare these equations with the usual Helmholtz dynamics. The variation of 
the Hamiltonian reads 

SH = f d 3 rv a (r)5v a (r) (121) 



Substituting here the velocity variation ( |113| ) we could drop the d a 8h term as it vanishes 
after integration by parts. As a result we find 

e aPl vp(X(p))^(p) (122) 



SX a (p) 

Finally, in the equation of motion we have 



d t X a (p) = - I d 3 p'5 3 (p - p'^y^Le^v, (X(p)) n v {p) (123) 



a 



Jaj3 



So, the Poisson brackets correspond to the motion in transverse direction to the gauge 
transformations. Or, to put it in different terms, we could modify the Helmholtz dynamics 
by adding the time dependent gauge transformations Sp(t) to the time shift of the vortex 
sheet 

d t X a (p) = v a (X(p)) + d a X a (p)d t 5p a (124) 

where 

n^ngMaHM) (125) 

Should we insist on the unmodified Helmholtz dynamics, we would have to admit that 
this cannot be achieved by any Poisson brackets. This is so called generalized Hamiltonian 



dynamics, where the formula (|122|) for the variation of the Hamiltonian cannot be uniquely 
solved for the velocity. The terms corresponding to the gauge transformation remain 
unspecified. In our opinion, this difference is immaterial, as the Helmholtz dynamics is 
indistinguishable from the one with the Poisson brackets. 
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